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Abstract 

We give some results about the dynamics of a particle moving in Euclidean three-space 
under the influence of the gravitational force induced by a fixed homogeneous circle. Our 
main results concern (f ) singularities and (2) the dynamics in the plane that contains the 
circle. The study presented here is purely analytic. 

In this paper we present a study about the movement of a particle in Euchdean three-space 
IR'^ on which the only acting force is the gravitational force induced by a fixed homogeneous 
circle. Even though this problem seems quite natural, we could not find in the literature any 
reference concerning the dynamics of it. All we could find were a few different ways of express- 
ing the potential function. Essentially all these expressions had already appeared in Poincare's 
Theorie du Potentiel Newtonien published first in 1899. Probably the most well-known ex- 
pressions are the one expressed in terms of elliptic integrals of the first kind and the one using 
the arithmetic-geometric mean given by Gauss. We believe this is the first study presented 
about the dynamics of this problem. 

Before we present our first result we need some notation and a couple of definitions. We 
want to study the movement in IR^ of a particle P under the influence of the gravitational force 
induced by a fixed homogeneous circle C. Denote by r = (x, y, z) G IR^ — C the position of the 
particle P. Also r = (x, y, i) denotes its velocity. According to Newton's Law the movement of 
P obeys the following second order differential equation: 

r= -VV^(r) (0.1) 

where V denotes the potential energy induced by C. The expression of V is given by V{v) = 
~ -fc \\r-ul ' '^'^^re A is the constant mass density of the circle C. We say that a solution r(t), 
defined in the maximal interval (a, b), has a singularity at b (or at a), if & < +00 (if a > —00). 
For p G ]R^, let dist{p,C) denote the Euclidean distance, inf^^cWP ~ ^11 1 from p to C. We can 
now ask: if lim^^f^- dist{r{t) , C) = 0, b < +00, does r(i) approach a well defined point in C ? A 
priori, r(i) could approach the circle without getting closer to any specific point in the circle. A 
singularity at c is called a collision singularity if there exists r* £ C such that r(t) r* , when 
t ^ c. Otherwise, the singularity is called a singularity without collision. Here is our first result. 



*The first author was partially supported by a CNPq/FACEPE research grant. The third author was partially 
supported by a research grant from CNPq, Brazil. 



1 



Theorem A. All singularities in the fixed homogeneous circle problem are collision singularities. 



In what follows we consider the fixed homogeneous circle C contained in the xy-plane and 
centered at the origin. Also, by rescaling we can consider the circle with radius equal to one 
(see section 1.3). Then the mass M of C is given by M = 2ttX. A quick examination of the 
problem shows that (see section 1.2 for more details) the z-axis, the horizontal plane (i.e. the 
xy-plane, which contains the circle) and any vertical plane (i.e any plane that contains the z- 
axis) are invariant subspaces of our problem. Also, any radial line (i.e. any line in the xy-plane 
that passes through the origin) is an invariant subspace. Here by an invariant subspace A we 
mean that any solution that begins tangentially in A is totally contained in A. The restriction 
of the problem to the z-axis is a one dimensional problem which is not difficult to treat. In 
fact it is a special case of a problem studied by Sitnikov (see [2], 0). Our two next results 
discuss the dynamics restricted to the horizontal plane. This problem is a central force problem. 
For r G IR — {1} define the function V{r) = V{ru), where u is any vector in the xy-plane of 
length one (note that we are using the same letter V for two different functions). The dynamics 
restricted to the horizontal plane has two cases: inside the circle and outside the circle. It is a 
classical result that using the polar coordinates (r, 9) of a point (x, y) in the horizontal plane it 

can be proved that svstem lHTTl in this case is equivalent to: 'r= "^~^^('")' ~5 (^^^ IS)- Here 
the constant K is the angular momentum of the solution. Hence this problem can be reduced in 
a canonical way to a problem with one degree of freedom (the variable r{t)) and once we know 
r(t), we can obtain 9{t) by integration. Note that K = ii and only if the solution lies in a 

d 

radial line (i.e. 6{t) is constant). Equation r= — — "^^^"^^ equivalent to 

r{t)^ -±U{T{t)), (0.2) 

where U denotes the effective potential U{r) = — + V(r). Note that U{r) depends on K, that 
is for each K we may have different functions U{r) (maybe we should write Uk{t) but we do 
not want to complicate our notation). Hence to determine the dynamics we have to know the 
behavior of the effective potential for all values of K. Note that outside the circle the effective 
potential is defined for 1 < r < oo and inside the circle the effective potential is defined for 
< r < 1 (for -fC / 0) and for -1 < r < 1 (for K = 0). 

Theorem B. The effective potential U (r) inside the circle has the following properties: 

i. -§pV(r) < 0, for < r < 1. Hence ^U{r) < 0, for < r < 1. 

ii. lim^._^i-U{r) = — oo. 

iii. For K ^ we have limj,^o+ U{r) = oo. 

iv. For K = we have that U{r), —l<r<l,is even. 

The following are direct consequences of the Theorem above. 

1. The graphs of U (r) for the cases K = Q and K ^ Q have the following form. 
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Fi gur6 0.1; Graph of the effective potential U inside the circle. 

Note that for K = 0, U{r), —1 < r < 1, has an absolute maximum value Eq = U{0) = V{0) = 
—M = —2itX , at r = 0. Hence r = or, equivalently r{t) = 0, t £ IR, is an equilibrium solution. 

2. By (i) of the Theorem, the force inside the circle is repulsive. 

3. As mentioned before, for K = a solution r(t) (with polar coordinates {r{t),6{t)) stays in a 
radial line, i.e. 6 is constant. Without loss of generality we can assume that this radial line is 
the X-axis, that is 6 = 0. Hence r{t) = (x(t),0, 0), where x{t) = r(t). Let E = E{x,x) denote 
the energy of r{t) = 0, 0). It follows that the phase portrait for = (i.e the level curves 
oi E = E{x, x) in the xx-plane) has the following form: 



Eo 




Figure 0.2: Phase portrait for K = 0. 

Consequently, depending on the energy the particle behaves in the following way. li E < Eq 
the particle comes from the point p = (±1,0,0) in the circle, stops before reaching the origin 
and then turns back to the point p. If E > Eq the particle comes from the point p in the circle, 
and goes all the way to the point —p. 





E<Eo 



E>Eo 



HE = Eq we have three cases: the particle comes from the point p in the circle and converges 
to the origin (in infinite time) , or comes from the origin (in infinite time) and converges to p, or 
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stays at rest in the origin. 




E = Eo E = Eo E = Eo 



4. For 7^ we have that a solution r(t) never passes through the origin. By Theorem B, this 
solution comes from the circle and returns to it (in finite time, by Theorem A). Note that there 
is a unique to such that r{to) is the closest point to the origin. Since r(to) = we have that 
r(to) is perpendicular to r(to)- We say that r(t) is normalized if r(to) hes in the positive y-axis. 
Note that, by symmetry, any solution (inside the circle) with K ^ can be obtained from a 
normalized one by a rotation. (See Fig. 0.3 below.) 

The following two Propositions complement the results given by Theorem B, for the case of 
a solution inside the circle. 

Proposition 1. Let r{t) be a solution of \U.ll contained in the horizontal plane and inside the 
circle. Let (a, b) be its maximal interval of definition. If r{t) approaches the circle as t ^ b~ or 
t — > a"*", then 

(1) b < oo or —oo < a, respectively. Hence, by Theorem A, the solution collides to a point in 
the circle, 

(2) when r{t) converges to the circle the speed converges to infinity and the velocity r{t) becomes 
orthogonal to the circle. 

It follows from this Proposition that the intervals of definition (a, b) of solutions inside the 
circle are as follows: (1) for K ^ we have that a and b are both finite (2) for K = and 
E ^ Eq we also have a and b are finite (3) for K = and E = Eq either a = — oo and 6 = oo or 
exactly one of a and b is finite. 

Proposition 2. Let r(t) be a solution of \().l[ contained in the horizontal plane, inside the circle 
and with angular momentum K ^0. Ifr{t) is normalized then the trace of the curve r(t) is the 
graph of an even convex function y = y{x). 
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Fi gurc 0.3; Trace of the normalized solution r(£). 
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We now state the results concerning the effective potential C/(r) for the case r > 1, that is 
outside the circle. 

Theorem C. The effective potential U{r), r > 1 (outside the circle) has the following properties: 

i. I: V{r) > 0. Hence ^{U{r) - |g) > 0. 

ii. limj.^i+ U{r) = —oo. 

iii. lirrir^oo U{r) = 0. 

Also there are Kq > 0, tq, 1 < tq < 2, and continuous functions ri = ri{\K\), r2 = r2{\K\) 
defined for \K\ > Kq, with 1 < ri < tq < r2 such that: 

iv. For K G {—Kq,Kq), U{r) does not have critical points. Hence from ii. ^ U{r) > and 
U{r) is an increasing function. 

V. For K ^ {—Kq,Kq) the critical points ofU{r) are exactly ri = ri{\K\), r2 = r2{\K\). 

vi. ^ U{r) > for r < ri and r > r2. Also ■^U{r) < for ri < r < r2. 

vii. ri{Ko) = r2{Ko) = tq. 

viii. Ti is d dGCTGdsifig fuTictioTi dfid Uttij^^^qq ri = 1. 
ix. r2 is an increasing function and limK^oor2 = oo. 

X. limK^ooU{ri{K)) = oo and limK^ooU{r2{K)) =0. 

Remark. The special value tq does not depend on the density A (or on the mass M) of the 
circle (see the proof of Theorem C). Since Kq = rQ-^V{rQ), the value Kq does depend on A. It 
is interesting to compare this with Corollary 11.81 

It follows from Theorem C that the graphs of U{r) have the following form: 
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Fi gure 0.4: Graph of the effective potential U outside the circle. 
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Note that for K ^ [— -K^O) Kq] the function U{r) has a local maximum at ri and a local minimum 
at r2- By (x.) of Theorem C, the local maximum value U{ri) at ri is eventually (as K ^ oo 
and ri ^ 1~) a global maximum and U{ri) +00. The following are direct consequences of 
Theorem C. 

1. From (i) of the Theorem, the force outside the circle is attractive. 

2. We discuss the stability of circular solutions. Here by stability we mean the following. A 
circular solution r(i) in the horizontal plane is stable if its reduced solution r{t) = ||r(t)|| = 
constant is a stable equilibrium position of the reduced problem r{t) = —■^U{r{t)). It is easy to 
see that a circular solution r(t) is stable iff the trace of solutions with initial conditions close to 
the initial conditions of r(t) stay close to the trace of r(t). Recall that for a fixed K, a circular 
solution with radius r has momentum K if and only if r is a critical point of U{r). Hence for 
K G (—KqjKq) there are no circular solutions with momentum K. For every K ^ [—Kq,Kq] 
there are exactly two circular solutions with radii ri(|Er|) and r2(|i^|). The circular solution with 
radius ri is not stable and there are "spiral" solutions approaching it from the inside and the 
outside. The circular solution with radius r2 is stable and has no "spiral" solutions approaching 
it. For K = Kq there is exactly one circular solution and it has radius tq. This circular solution 
is unique among all circular solutions. It is not stable and have spiral solutions approaching it 
only from the inside. Also it can be approximated by stable circular orbits from the outside 
(with varying momenta K). It follows that a circular solution with radius r is stable if and only 
is r > rg. Note that a non-circular solution with momentum Kq either collides to the circle, 
escapes to infinity or is the solution that approaches the circle (in infinite time). 

3. Let E = E{r,r) denote the energy of r{t) = ||r(i)||. Using Theorem C we can deduce the 
phase portraits. We have three cases: K G {—Kq,Kq), K = Kq, K ^ [—Kq,Kq]. For the first 
two cases the phase portraits are shown in figures 0.5 and 0.6 respectively, and the solutions 
have the following behavior: 

(i) K £ {—Kq, Kq). In this case all solutions either collide to the circle or escape to infinity. 

(ii) K = Kq. In this case all solutions (different from the unique circular solution of radius tq) 
either collide to the circle or escape to infinity or converge to the circular solution with radius 
vq (in infinite time). 



E<U(r„) U(ra)<E<0 




E<0 



E = U(r,) 



Figure 0.5: K e {~Ko,Kq). 



Figure 0.6: K = Kq 
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(iii) We analise now the case K ^ [—Kq,Kq\. Let E = U{ri). We have also three cases: 
E > 0,E < 0,E = (see figures 0.7, 0.8, 0.9). Note that, by (m.), (mi.) and (x.) of Theorem 
C all three cases do happen. For E' > we have (see figure 0.7): In region I the solutions 
are bounded, and converge to the circle. In region II the solutions are unbounded, come from 
infinity and converge to the circle. In region III the solutions are unbounded, come from the 
circle and escape to infinity. In region IV the solutions are unbounded, approach the circular 
solution (with radius ri) and escape to infinity. In region V the solutions are bounded, and stay 
close to the circular solution (with radius r2). 

In figure 0.7 the solution represented by 71 comes from the circle and converges to the 
unstable circular solution (with radius ri). The solution represented by 72 comes from the 
circular solution (with radius ri) and converges to the circle. The solution represented by 73 
comes from the circular solution (with radius ri) and escapes to infinity. The solution represented 
by 74 comes from infinity and converges to the circular solution (with radius ri). 

For E = we have also three cases: E < E, E = E, E > E. We have the following analysis 
(see figure 0.8 below): For E < E in region I the solutions are bounded, and converge to the 
circle; in region IV the solutions are bounded and stay close to the circular solution (with radius 
r2). For E > E we have: In region II the solutions are unbounded, come from the circle and 
escape to infinity. In region III the solutions are unbounded, come from infinity and converge 
to the circle. For E = E we have: The solution represented by 71 comes from the circle and 
converges to the unstable circular solution (with radius ri). The solution represented by 72 
comes from the circular solution (with radius ri) and converges to the circle. The solution 
represented by 73 comes from the circular solution (with radius ri) and escapes to infinity. The 
solution represented by 74 comes from infinity and converges to the circular solution (with radius 

n). 

For E < we have four cases: E<E,E>0, E<E<0, E = E (see figure 0.9 below). 
For the case E < E in region I the solutions are bounded and converge to the circle; in region 
IV the solutions are bounded and stay close to the circular solution (with radius r2). For E > 
we have: In region II the solutions are unbounded, come from the circle and escape to infinity. 
In region III the solutions are unbounded, come from infinity and converge to the circle. For 
E < E < the solutions are bounded, come from the circle, approaching the circular solutions 
and return converging to the circle. Finally, iov E = E we have: The solution represented by 
7i comes from the circle and converges to the unstable circular solution (with radius ri). The 
solution represented by 72 comes from the circular solution (with radius ri) and converges to the 
circle. The solution represented by 73 comes from the circular solution with radius ri, approach 
the circular solution with radius r2, and return converging to the circular solution with radius 

n- 




Figure 0.7: £; = c/(ri) > 0. Figure 0.8: £; = c/(ri) = 0. Figure 0.9: s = [/(n) < 0. 
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We mentioned above that a circular solution of radius r in the horizontal plane is stable 
(among solutions in the horizontal plane) if and only if r > tq. The following question arises: 
are these solutions stable among all solutions in IR^? The answer is affirmative, but before 
we state this result we have some comments. Since our original problem (in IR^) is invariant 
by rotations around the z-axis, the problem can be reduced in a canonical way to a problem 
with two degrees of freedom (see section 5). In fact, every solution r(t) = (x{t),y{t), z{t)) can 
be written as r(t) = (r{t)cosip{t),r{t)si7up{t), z{t)), i.e the cylindrical coordinates of r{t) are 
(r(t), (^(t), z(t)), and {r(t),ip{t), z{t)) satisfy a system of equations (system 1)5. 7p ). Moreover, 
once we know r(t) we can obtain ip{t) by integration. We say that (r{t),z{t)) is the canonical 
projection solution of r(i). That is, the canonical projection solution is obtained from the solu- 
tion by the map {x,y,z) i— > (r, z), r = \/x'^ + y^. We have that {r{t),z{t)) satisfy the first two 
equations of (|5.7|) . Note that under this projection the circle C projects to a point xq with rz 
coordinates (1,0). We extend the definition of a circular solution (for central force problems) to 
this context in the following way: We say that a solution is circular if its canonical projection 
solution (r{t),z{t)) is an equilibrium position of the system formed by the first two equations 
of (|5.7|) . Also, we say that a circular solution is stable if its canonical projection solution is a 
stable equilibrium position of the system formed by the first two equations of ()5.7|) . It is easy 
to see that a circular solution r{t) is stable iff the trace of solutions with initial conditions close 
to the initial conditions of r{t) stay close to the trace of r(t). 

Proposition 3. All circular solutions lie in the horizontal plane. Moreover, a circular solution 
(in the horizontal plane) of radius r is stable in IR^ if and only if r > tq, where tq is as in 
Theorem C. 

We present two more results. Fix A > 0. For e > denote by W{r, e) the potential function 
induced by a the fixed homogeneous circle in the xy-plane centered at (^,0,0) with radius ^ 
and density A. We write W{x, z; e) for the restriction of this potential to the xz-plane. Define 
VW{x,z;0) := 64 A^^, {x,z) / (0,0). Let VW{x,z;e) be the gradient of W{x,z;e) with 
respect to the variables x and z. Also, let A = {{x,z;e); (0,0) 7^ {x,y) / (|,0)}. Thus 
VW{x, z; e) is defined on A. 

Proposition 4. VW is continuous on A. 

Note that VW{x, z; 0) = GiXVlniVx"^ + z^) and recall that 2Xln{Vx^^Tz^) = \ln{x'^ + z^) 
is the potential induced by the infinite wire (with constant density A and infinite mass) orthogo- 
nal to xz-plane intersecting the xz-plane at the origin. Hence by Proposition 4 the problem of the 
fixed homogeneous circle with large radius, constant density A, and centered at (i, 0, 0)) can be 
regarded as a perturbation of the problem of the infinite homogeneous straight wire with density 
32A. (To obtain the same density A, instead of 32A, it is enough to take circles with radius ^.) 
This result is a key element in the proof of the existence of periodic orbits near the circle (see |IJ). 

There are a few textbooks where it is claimed that the potential function of fixed homoge- 
neous circles with large radii "converge" to the potential function of the infinite straight wire, 
but we have not found satisfactory proofs of these claims. In fact, it seems that in order to avoid 
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the problems with the infinities, we have to work with the gradients, instead of the potentials. 
This necessity seems to be related to the fact that if we try to calculate the potential of the 
infinite straight wire by integrating \\p-(x,o o)\\ obtain infinity, but the integral for the 
gradient does converge (or can be computed easily from Gauss formula). 

Finally we present our last result. The potential V of our problem at P can be written 
as V{P) = —^^^^5 where a{D,d) is the arithmetic-geometric mean of D > and d > (see 
Remark 1.1). Here d, D are the maximum and minimum distances D = D(P), d = d{P) of the 
point P to the circle. As we mentioned above, this formula of V{P) was given by Gauss. In our 
last result we give a formula for {-^V, ^^) in terms of D,d and their successive geometric and 
arithmetic means dj , Dj . 



Proposition 5. We have the following formulas: 



where x^j:Sii,ilg-nU{l 



This paper has 7 sections. In section 1 we give some preliminary results that will be needed 
later. These include some basic facts about symmetry, elementary properties of the potential 
function as well as some identities related to the potential function. In sections 2, 3 and 5 we 
prove Theorems A, B and C, respectively. At the end of section 5 we prove Proposition 3 and 
also prove that the only equilibrium position is the origin. By symmetry it is clear that the 
origin is an equilibrium position but it is not too trivial to prove that it is the only one. In 
section 4 we prove Propositions 1 and 2. Finally in sections 6 and 7 we prove Propositions 4 
and 5, respectively. 



Preliminaries. 



1.1 The Potential. 



As before consider the fixed homogeneous circle C with constant 
density A and radius p, contained in the xy-plane, and centered at 
the origin. The mass M of C is 27rXp. We denote by P = (x, y, z) 
the coordinates of the position of a particle in IR^ \ C and let V{P) 
be the potential at P induced by C. Let D = D{P) and d = d{P) 
denote the maximum and minimum distances from the particle P 
to the circle. We have the following expressions for these distances: 



have that (see |7|, p. 196) 



y2+p)2 + 2;2 



and d'^ 



(v^c2 + y2 . 



Also we 



V{P) = -4Xp 



d^j 




\/ d'^cos'^ip + D'^si'nP'ip 



Figure 1.10: Fixed circle cen- 
tered at the origin. 

(1.3) 



9 



Define T : {(x,y); x > 0, y > 0} ^ IR by T{d,D) = /^^ ."Ir^.-.r Then we have 

V{P) = —4:\pT{d, D). By the change of variable ip = ^ — 6, it can be proved that T is 
symmetric with respect to d and D, that is, T{d,D) = T{D,d). It was proved by Gauss (see 
0, p. 197-199) that V{P) = —-^^jyj where a{d,D) is the arithmetic-geometric mean of d and 
D. Rence T{d,D) = ^^. 

Remark 1.1 Recall that the arithmetic-geometric mean a{m, n) of two positive numbers m, n is 
defined in the following way. Set mi = and ni = ^Jmn and define inductively the sequences 
{mj}, {nj} by mj+i = HH+Hi^ j^j+i = ^miTii. We have m > mi > m2 > ... > n2 > ni > n. It 
can be shown that lim rrij = limnj. This common limit is called the arithmetic- geometric mean 
a{m, n) of m and n. Note that ifm = n, then a{m, m) = m. Note also that a{m, n) = a{mi,ni). 
Since T{m,n) = 2a-(m n) ' have T{m,n) = T{mi,ni). 

Define now f{t) = , J\ . < t < 1. Then T{d,D) = = ^ ,"1 < 

d<D. Hence V{P) = -^f (-g^) . It follows that / = 2^^, 0<d<D. Setting L> = 1 
we have that f{d^) = 2a(d,i) ' ^ — 1- 

Remarks. 



the first kind. 



(1) f(t) = / ' , . = / ' ^ = = K(VT^), where K is the elliptic integral of 

~ ^)"^ • This series can be obtained from the power series 




(2) fit) = ^ . ^ ^ 

\ n=l 

oiK (see 0, p. 197). 

Lemma 1.2 (i) f is a decreasing function and f is an increasing function, 
(ii) a{d, l)>Vd 
fit) < f ^ 

H \f'{t)\ < W 

(v) \f'{t)\ < f ^ 

(vi) limi^o+ fit) = +0O 

(vii) limt_o+ tf'{t) = -\ 

Proof, (i) f'{t) = -\ f^ (cos^7+t!i%e f'^ ^ '^^ich imphes that / is a decreasing function. 
fit) = 4 - — ^'"''^'^^ > 0, which implies that f is an increasing function. Note that 

/(")(t) < 0, if n is odd, and /(")(t) > 0, if n is even, [ii) By the definition of a, we have 
that a{d, 1) > (geometric mean of d and 1) = ^/d. {Hi) Setting D = 1, d^ = t, and applying 
(ii), we obtain (Hi), (iv) i , ,f "''^'^L.a,, = ^ \^ \-, — tt^-^t^i — iwrr^dO < 

i^^, (because (^cos'^e+tsin'^e) < ^, t > 0.) The proof of (u) follows directly from (Hi) and 

(iu). (ui) Since VcosW+tsEFe < VcosWTt < cose + Vt we have f{t) > /g' ^^J^^ . A 



10 



change of variable shows /g' 

TT TV 



de 



k 







^ fa 

^ + Vi)-ln{Vi) 



—z — 7=. Hence, since sinO < 9 we have f(t) > 
^ cxD as t ^ 0"^. (vii) We have the following power 



series for the elliptic integral of the first kind K{k) (see [7|, p. 203), K{k) = 



1 + f + 64 

taking k\ 



9 1,4 



In-, 



4 64 ^ ■ • • 
-1.- 



t we have f{t) 
1 



4 128 
^ 4 ^ 64 ^ • • • 



, where /cf 



4 ^ 128'' ^ 



+ 



1 _|_ 1 _|_ -9-/-2 I 
^ 4 ^ 64 '' ^ 



Since /(t) = K{^/T^) 
... It follows that /'(t) = 
therefore limt^of' {t)t = 



L _i_ _i_ 
4 128'' 



1.2 Symmetries of the potential. 

In this section we study the symmetries of the fixed homogeneous circle problem and determine 
the invariant subspaces. As before, we are considering the fixed homogeneous circle C, contained 
in the xy-plane and centered at the origin. We call the xy-plane (which contains C) by horizontal 
plane, and any plane that contains the z-axis by vertical plane. 

A symmetry ip of our problem is a transformation ij) : IR^ — C ^ IR^ — C that preserve V, 
i.e V is invariant by ip, that is = V. The symmetries of V are evident from the geometry of 
the problem: they are the rigid transformations of IR^ that preserve C: 

Proposition 1.3 The potential V is invariant by: 
(i) Rotations about the z-axis, 
(a) Reflection with respect to the horizontal plane, 
(Hi) Reflections with respect to a vertical plane. 

Proof. It follows from (|1.3j) and from the expressions of D and d (see section 1.1). ■ 

Corollary 1.4 The gradient field is invariant by: 
(i) Rotations about the z-axis, 
(a) Reflection with respect to the horizontal plane, 
(Hi) Reflections with respect to a vertical plane. 

Proof. Since all symmetries above are isometrics of IR^ \ C (with the canonical fiat metric), 
the Corollary follows from the Proposition above and from a classical result in the elementary 
theory of mechanical systems on Riemannian manifolds. ■ 

Corollary 1.5 If v{t) is a solution of hU. then %l)r{t) is also a solution of hU.l\) . where is 
one of the symmetries above. ■ 

Corollary 1.6 The z-axis, the horizontal plane and the vertical planes are invariant subspaces, 
that is, a solution of \U.1\) that starts in one of these spaces, with velocity contained in it, stays 
there for all the time in which the solution is defined. 

Proof. It follows from the following facts: 
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- the z-axis is the invariant space of ah the symmetries in (i) of Corollary 11.41 

- the horizontal plane is the invariant space of the reflection in (ii) of Corollarv 11.41 

- a vertical plane is the invariant space of a refection in (Hi) of Corollarv 11.41 ■ 

Note that a radial line (i.e. any line in the xy-plane that passes through the origin) is also 
an invariant subspace since it is the intersection of the horizontal plane with a vertical plane. 



1.3 Properties. 

We will need to consider circles with variable radius and mass. Write V{r,p,M) to denote the 
potential induced by C, contained in the xy-plane and centered at the origin, with radius p and 
mass M, and Vy(r, M) to denote the gradient (with respect to r) of V{r, p,M). 

Lemma 1.7 The potential V of the fixed homogeneous circle problem satisfies the following 
identities: 

(i) V{r,p,cM) = cV{r,p,M), for c e IR, 

(ii) V{cr,cp,M) = ^V{r,p,M), for c> 0, 
(Hi) W{cr,cp,M) = j^W{r, p,M), for c>0, 
(vi) VV{r,p,cM) = cVV{r,p,M), for c £ IR. 

Proof. It follows directly from the definition of V{r,p,M) = — ^ Jq^ — where e*^ = 
{cos6, sin9, 0). ■ 

Corollary 1.8 Let p,(^,M,N be positive numbers. Ifj'{t) is a solution of r (t) = — Vy(r,/9, M) 
then s{t) = ^r [^Jj^ is a solution of s {t) = -VV{s, (, N). 

Proof. It follows from Lemma 1 1.71 by a direct calculation. ■ 

Remark 1.9 Note that if r{t) and s(t) are as above, then they have the same qualitative prop- 
erties. For instance, if r{t) is periodic then s{t) is also periodic (with period -^==^, where T is 
the period of r(t)). 

Note that Lemma 11.71 and Corollary 11.81 imply that in the study the fixed homogeneous 
circle problem we can assume the mass and the radius to be equal to one. 



2 Singularities: Proof of Theorem A. 

To prove Theorem A we first show the following Proposition, which is of general nature. In the 
next Proposition dist denotes "Euclidean distance" . 

Proposition 2.1 Let F : 17 — > iR, C IR" open. Let r{t), t S (a, 6), be a solution of r = 
—W{r) such that there exist: 

(1) vo,vi e M, vo < vi, with dist{V~^{vo),V~'^{vi)) > 0, 
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(2) ti < Si < t2 < S2 ■ ■ ■ , ti, Si £ {a,b), i £ IN, satisfying the following two properties: (a) 
V{v{ti)) = vq, and (h) V(r(sj)) = vi. Then b = +00. 

Proof. Since V is continuous, for every interval [ti,Si] we can choose, for all i, an interval 
[ths*] ^ [ti,Si] such that vo = V{r{t*)) < V{r{t)) < V{r{s*)) = vi, for all t e [t*,s*]. 

Since the total energy E = E{r{t)) = V{r{t)) + ^||r(t)p is constant, we have that ||r(t)|| = 
y/2{E - V), with V = V{r{t)) and E constant. Then we have ||r(t)|| < ^y2{E - vq), for t E 

s*]. Since the length of the curve r{t) between [t*, s*] is larger than the distance 
d = dist(y~^{vo), V~^{vi)) > 0, we have 

d < jl^ \\r{t)\\dt < p ^2{E-vo)dt = {s* - t*)^2{E-vo) < is^ - ti)^2{E-vo). 
Hence, . < (sj — ti), for all i. Since we have an infinite number of disjoint intervals 

\/2(E-vo) 

(ti, Sj), with < , < (sj — tj) for all i, we conclude that b = +00. ■ 

In what follows of this section let C be the fixed homogeneous circle in IR^ centered at the 
origin, contained in the horizontal plane and with constant density A and radius 1. Also let V 
be the potential induced by C and r{t) be a solution of !•= — Vy(r), defined in the maximal 
interval {a,b). 

Proposition 2.2 If b < +00 then limi^f)-dist{r{t),C) = 0. Analogously, if a > —00 then 
limi_,a+dist{r{t),C) = 0. 

To prove this Proposition we need the following Lemmas. Let {rn}„g]N C IR^ — C. 
Lemma 2.3 /im„_^+ooy(r„) = if and only if linin-^^ooW^nW = +00. 

Proof. Suppose first that ||r„|| — > +00. For all u £ C, we have = p. Hence, for ||r|| > p, and 
ti G C, we have ||r — u\\ > |||r|| — ||n||| = |||r|| — p\ > ||r|| — p. Hence < —V{r) = \\r-u\ — 
Jq ji-i_p du = p^jz^) where M = f^Xdu is the mass of C. It follows that /im„_^+oo^(rn) = 0. 

Suppose now that V{vn) 0. Since ||ii|| = p, for u £ C, we have j^Y-u\ — \r\+p • Conse- 
quently, < = jrf+^du < Jc p^d^i = -V(r). Therefore, /im„_^+oo||r„|| = +00. This 
proves the Lemma. ■ 

Lemma 2.4 limn^+ooV {rn) = —00 if and only if limn^+oodist{rn,C) = 0. 

Proof. Suppose that limn^+ocy{T^n) = —00. Since < dist{rn,C) < ||r„ — for all n € C, we 
have < -y(r„) = pH^dn < oHt^du = M^^^^^y It follows that dist{Tn,C) ^ 0. 
Conversely, suppose that dist{rn,C) — > 0. Let d, D be as in section 1.1. We have d(r„) = 

dist{rn,C) and L>(r„) 2p. It follows from (m) of LemmaOthat y(r„) = fiw^) 
—> —00. ■ 

Proof of Proposition 12.21 First note that, by the two Lemmas above, y~^(c) is compact, 
for all c < 0. Now, let r(t), t £ {a,b), b < +00, be a maximal solution of !■= — Vy(r). 
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Then limf^i,-V {r{t)) exists (it could be finite or infinite), otherwise, by the continuity of V we 
could choose two sequences (tn), (sn) such that ti < si < t2 < S2 < ■ ■ ■ , U, Si G {a,b), with 
V{r{ti)) = vo,V{r{si)) = vi, vq < vi. Since V^^{vq) and V^^{vi) are disjoint compact and 
non-empty, we have that distCV'^^ {vo),V^^ (vi)) > 0. Hence, by Proposition 12. II we would have 
that b = +CXD, a contradiction. 

Since imageV C (— oo,0), it follows that limi_,ij-V{r{t)) is either 0, or it is a number 
V* 7^ 0, or it is equal to -co. We show that the two first possibilities do not happen. It 
will then follow that linii^f,-V{r{t)) = —oo and therefore, by Lemma 12.41 we would have that 
limt^i,-dist{r{t),C) = 0, which proves the Proposition. 

If /imj_^j- y(r(t)) = 0, then there exists Iq such that for all t > Iq, vq < V{r{t)) < 0, with 
Vq = V'(r(to))- On the other hand by Lemma 12.31 we have that ||r(t)|| — > +oo when t — > b~ . 
Moreover, since E = V{r{t)) + |||r(t)|p, E is constant then ||r(t)|| < y/2{E - vq), for all t e 
(to,b), that is, the velocity is bounded in this interval. We also have that J^* ||r(t)||(it > \\r{t) — 

r{to)\\ > dist(r{t),V-^{vo)) =: df Hence, dt < {t - to)^2{E-vo), and then, to + -j^ = < 

t, for ah t £ {to,b). 

Since V~^{vo) is compact, there exists r > such that V~^{vo) C B{0, r) and since ||r(t)|| — > 
+00, given n > 0, there exists such that for all t G [tn,b), r{t) ^ B{0,n). This implies that, 
dt = dist{r{t), V~^{vq)) > n — r, for all & > t > t„. In particular, dt„ > n — r, and for b > tn > to, 

we have b > tn> tQ-\ — . """^ Hence, lirrin-^oo I H — , ""^ I = +cxd, and we conclude that 

^^2{E-vo) \ y/2{E-vo) J 

b = +00, a contradiction. 

Finally suppose that limt^i,-V{r{t)) = v*, with v* £ {—oo,0). Hence, given e > 0, there 
exists te such that for all t > t^, V{r{t)) £ (v* — e,v* + s). Let vq = v* — e, vi = v* + e; we 
can suppose that < e < \v*\, hence vi < 0. Clearly the set l^~^([fo, fi]) is not empty and 
since V is continuous, this set is closed in IR^ \ C. Moreover, by Lemma 12.41 this set is closed 
in IR^, and by Lemma it is bounded (because vi < 0). Hence, ^"^([t'o, wi]) is compact and 
r{t) e V-'^{[vo,vi]), for ah t > te- 

On the other hand, from the energy equation, we have ||r(t)p = 2{E — V{r{t))), therefore, 
for all t > t^, we have ||r(t)|| < \/2{E — fo) = C2. Hence the maximal solution (r(t),r(t)), 
t G (a, of the system of first order differential equations 

J r = V 

I V = -vy(r) 



is contained in the compact y~^([T;o, I'l]) x -6(0, C2), for t G {t^,b). It follows from a classical 
result in the elementary theory of differential equations (see jH]) that b = +00, a contradiction. 
This proves the Proposition. ■ 

Proof of Theorem A. Let r(t), t € (a, b), b < +00, be a maximal solution of r = — Vl^(r). By 
Proposition 12.21 we have that lirrit^i,- dist{r{t) ,C) = 0. To prove that limt^i,-r{t) = r* G C, we 
write this system in cylindrical coordinates (r, 9, z) in (IR^ \ {^ ~ axis}). We have (see equations 
1)5. 7() of section 5): 

.. _ K'^ dV .. _ dV 
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with 6 = , where K is constant and V(r,z) = V(rcos9, rsin6,z). 

Since the circle in cylindrical coordinates is given by C = {{p,(p,0),ip G IR}, we have that 
showing lim^^fj-r{t) = r* G C is equivalent to showing lim^^f,-r(t) = p, lim^^i,- z{t) = and 
limf^fj-6{t) = for some ^0- The two first limits follow from the fact that /imj_,b-(i«st(r(t), C) = 
0. We will now prove that lim^_^i^-9[t) = 6q. If = 0, 9{t) is constant, and we have nothing to 
prove. Suppose then that K > 0. Hence 9{t) is an increasing function. Thus, to prove that the 
limit of 9{t) exists, it is enough to prove that 9{t) is bounded above, for t in a neighborhood of 
b. 

/■* K 

Since r(t) p, there exists to such that, for t > to, r(t) > f . Hence 9(t) = / „ ds + 

Jto r{sy 

ft 4^ 4:K(b — tn) 

9(to) < / -^ds + 9(to) < — — + 9(to) < +00 for ah t G (to, b). Therefore, the limit of 

Jto _ P^ 
9{t) when t ^ b~ exists. This proves the Theorem. ■ 



3 The Dynamics Inside the Circle: Proof of Theorem B. 

In this section we consider again the fixed homogeneous circle C contained in the xy-plane, 
centered at the origin and with radius 1 and constant density A. Let r = {x, y) be the position 
of the particle in this plane, under the infiuence of the gravitational attraction induced by C. 
Let r = ||r|| and 9 be the polar coordinates of r = {x,y). Also, let D and d be as in section 
1.1. In the horizontal plane we have = (r + 1)^ and = {r — 1)^ = (1 — r)^. Hence the 
expression for the potential H1.3() becomes 



y(r) = V{r) = -4A / 

JO 



d9 



\/(r + lYcos'^9 + (r — lYsin^9 

Since V depends only on r, we have a central force problem and (see j3|, ch.3) the system 
of equations r = — Vy(r) (restricted to the horizontal plane) is equivalent to the system: 



(3.4) 

K 

where = (r x r, 63) = xij — xy \s the angular momentum. Note that r = ^ — ^^C'^) ~ 
d ( -K^ V{r) \ which is equivalent to 



dr \ 2r'^ 



-^Uir), (3.5) 
ar 



where U(r) = — ^ + V(r). The function U(r) is called the effective potential . 
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Proposition 3.1 For < r < 1, ^V^(r) < 0. 

Proof. In the horizontal plane and inside the circle we have D = r + 1 and d = 1 — r. Hence 
Di = 1 and di = \/l — r"^, where Di and di are the arithmetic and geometric means of D and 
d, respectively. Therefore, by Remark II. II of section 1.1, we have 

V{r) = V{r) = -4AT(1 + r, 1 - r) = -4AT(1, Vl-r^) = -4A H -=J^==. (3.6) 

Jo V 1 — r^sm^d 

TT -On 

Differentiating we have j-V^(r) = -4A / ^ dO < 0. ■ 

Proof of Theorem B. Differentiating we have -^U{r) = — ^ + -^V{r). Hence, by the 
Proposition above we have ^U{r) < for < r < 1. This proves (i) of Theorem B. 

When r ^ 1~ we have that the particle tends to the circle. Hence, by Lemma 12.41 V{r) 

tends to — oo. Consequently U{r) = ^ + V{r) tends also to — oo. This proves (ii). Note that 
{Hi) and {iv) follow from the definition of U{r) and 1)3. 6|) . ■ 



4 Proofs of Propositions 1 and 2. 

First we prove Proposition 2, but before we give some remarks. Let r(t) be as in the statement 
of Proposition 2. Without loss of generality we can assume that to = (see definition of a 
normalized solution). 

(1) Recall that the curvature of the curve r(t) is by definition 
where n(t) = ^^^^ and r^(t) = ||r(t)||n(t) = (-y(t), x(i)). 

(2) Note that k{t) is continuous and is never zero because if k(t) = for some t, then i'(t)_Ln(t) 
and follows that r{t) is radial (because r{t) is always radial), and then the solution r(t) would 
be radial, that is -ftT = 0, a contradiction. 

(3) r(0) _Li:'(0) and i:'(0) points upward. It follows that n(0) = cr{0), for some c > 0. Then 
A;(0) = c||i:^(0)p > 0, with c > 0. In this case, since k is continuous and is always non-zero, we 
have that k > 0, for all t. 

Proof of Proposition 2. Let r{t) = {x{t),y{t)) be as in the statement of Proposition 2 and 
recall that we are assuming = 0. Let (a, b) be the maximal interval on which r(t) is defined. 
Then a < < 6. 
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Claim. x{t) / 0, for all t G (a, 6). 

We prove that x{t) / 0, for all t G [0,6), the proof for t G (a, 0] is similar. Suppose 
that there exists to G [0,6) such that x{to) = 0. Let to = min{t > 0; x{t) = 0}. Since x 
is continuous and x(0) > 0, it follows that to > 0. Then x{t) > 0, for t G [0,to)- We have 
r(to) = (x(to)) y(*o)) = (0,y(to))) with y(to) / 0, because r{t) ^ 0, for all t. We can write 
r(t) = a{t)(cosif{t), sinip{t)), witha(t) = ||r(t)|| > and (/?(0) = 0, continuous. Differentiating: 
r{t) = a{t){cos(p{t),sin(p{t)) + a{t)(p{t){—sin(p{t),cos(p{t)). Hence A; (t) = { r(t), ^'^{t) ) = 

° a^^(t{*^ = ^1^. Because the curvature is positive (see Remark (3) above) this shows that 
ip{t) > and ip{t) is an increasing function. Since x{to) = 0, ip{to) = f (hence y{to) > 0), 
or (f{to) = ^ (hence y{to) < 0). Since is an increasing function we have that (p{to) = ^, 
which implies that y{to) > 0. Since x is increasing on (0,to)) we have that x{to) > 0. Recall 
that r{t) is radial and expansive, that is, r{t) = b{t)r(t) = b{t){x{t),y{t)), with b{t) > 0. Hence 
k{to) = b{to) {{x{to),y{to)),{-y{to),0)) = b{to){-x{to)y{to)) < 0, a contradiction. Therefore 
there is no to such that x{to) = 0. This proves the claim. ■ 

It follows from the claim that x{t) > 0, a < t < b. Hence x{t) is a increasing function. In 

this way the function t x{t) is one-to-one and it follows that x{t) possesses an inverse t = t{x). 
Define f{x) = y(t{x)). Note that the graph of / is equal to the trace of r. Differentiating / 
with respect to x, we have 

Differentiating again we have 

^/(.) = ^(r,||r||n)=^>0 

because x > and A; > 0. Hence the trace of r(t) is given by the graph of a convex function 
/. By the symmetry of the problem, the solution rit) is symmetric with respect to the y-axis. 
Hence the function / is even. ■ 



Proof of Proposition 1. Let r{t) and (a, 6) be as in the statement of Proposition 1. Let r{t) = 
||r(t)||. Then r{t) satisfies r = —-^U{r). Suppose that r{t) —>■ 1" , when t —>■ b~ . By Theorem B 
we have limt^f,-U{r{t)) = limr_i-U{r) = — oo. Thus there exists to such that U{r{t)) < E — 1, 
for t G [to,b), where E = ^r^ + C/(r) is the energy of r{t). Hence, for t > to, r{t) ^ 0. Moreover, 
fit) > (because r(t) 1~). Consequently, r{t) is one-to-one on [to,b), and has an inverse 
t = t{r), To < r < 1, ro = r{to). Prom E = ^f'^ + U{r), we have that r = ^/2{E - U{r)). Then 

dt = /; , and it follows that tir) = T , '^^ +to, ro < r < L In this way 

6 = lim,^^-t{r) = l\ ^2{E-U{r)) + *o < -^-0 :S + *o - 72^1 - ^o) + *o < +°°- ™^ P™^^^ 
part (1) of Proposition 1. 

We now prove part (2) of Proposition 1. Again, let r(t) = {x{t),y{t)) and (a, 6) be as in 
the statement of Proposition 1. Assume ||r(t)|| —>■ 1, when t —>■ b~ . We have E = ^||rp -|- V{r). 
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Then ||rp = 2[E — V^(r)). When t ^ b~ , ||r|| — > 1 hence V{r) — oo. Therefore ||r|| +00. 
Without loss of generahty suppose that r{t) converges to the point (1,0) of the circle, that is 
x{t) — > 1 and y{t) — > when t ^ b~ . 

We shall show that r{t) becomes horizontal as t — > that is when t — > 6~. First, 

note that x{t) / for t close to b. Moreover, x{t) +00 when t ^ 6^. To see this suppose that 
there exists a sequence t„ — > b~ with |i;(tn)| < M for some M. Then \y(tn)\ +00 (because 
||r(^n)|| — *■ +00). Hence \K\ = \x{tn)y{tn) — y{tn)x{tn)\ — > +00, a contradiction because K is 
constant. 

Now, suppose that lirrii^^- |f|^ 7^ 0. This implies that there exists a sequence with 
tn — *■ b~ such that > 5, for some 5 > 0. Hence, we have \K\ = \x{tn)yitn) - y{tn)i{tn)\ > 

\x{'tn)\\y{tn)\ - \y{t„)\\x{tn)\ > \x{t,^)\{S\x(t„)\) - \y{t„)\\x{tn)\ = \x{tn)\ (^5\x{u)\ - |y(tn)|^. Talclug thc limit 
when tn — > we have that \K\ +00, a contradiction. Therefore, limt—,^ = 0. ■ 



5 Dynamics Outside the Circle: Proof of Theorem C. 

As in section 3 we consider the fixed homogeneous circle C contained in the xy-plane, centered at 
the origin and with radius 1 and constant density A. Also r will denote the position of a particle 
in this plane, under the influence of the gravitational attraction induced by C. Let r = ||r||. 

Proof of Theorem C. Let D = D{r) and d = d{r) be as in section 1.1, that is they are the 
maximum and minimum distances from r to the circle. In the horizontal plane and outside the 
circle we have D = r + 1 and d = r — 1. Hence Di = r and di = ^Jr^ — 1, where D\ and d\ 
are the arithmetic and geometric means of D and d, respectively. Therefore, by Remark 1 1.1 1 of 
section 1.1, we have 

Viv) = Vir) = -4Anr+l,r-l) = -4A T(r, J - 1) = -4A /' '^^ — = -4A 

io \^W^TF^I)We io Vr^ - sin' 

Differentiating we have 

^no-4Ar- — ^-^d^ > 0. 

This proves {%) of Theorem C. When r — > l"*" we have that the particle tends to the circle. 
Hence, by Lemma f2.41 V{r) tends to —00. Consequently \J(r) = — ^V{r) tends also to —00. 
This proves {ii). Note that {iii) follows from the definition of \J(r) and Lemma 12.31 

We now prove {iv)-{ix) of the statement of Theorem C. Since by definition U{r) = §7 + V{r) 
we have that ■§pU{r) = if and only if = ■^V{r), or equivalently = r^-^V{r). Define 

3/2 

- fr2 - sin'^e) 



18 



Lemma 5.1 The function g{r) has the following properties 

a. lim^._^i- g{r) = oo, 

b. linir—oo ^ = 27rA, 

c. g{r) has exactly one critical point tq, 1 < tq < 2. 

Proof. If r ^ 1~ the particle approaches the circle therefore, by Lemma 12.41 V{r) — oo. 
Hence -^V{r) is not bounded when r ^ 1~ . But -£7V{r) - 4A Jo j:^z^^^^0j2 dO < therefore 
^y(r) ^ oo, as r ^ This proves (o.). For (6.) we have ^ = 4Xj^ ^^2_",-„2g^3/2 dd = 
^0 (1 a.t^t))3/2 '^^ ^ ^'^i ^ ^ttA. To prove (c.) we compute the first two derivatives of g{r). A 
direct calculation shows: 



{r-2-si-n?e\ •'° ( 



It follows that ^ g(r) > 0. This together with (a.) and (b.) imply that g{r) has a 
unique minimum at some point rg > 1. Finally we show that tq < 2. For this just calculate 
^ g{2) = 4A8r/ ,/-^T'fM de = 32A ^''"'CL.^ dO > 0. ■ 



Note that, by definition of g, tq does not depend on A even though g(r) does. From the 
Lemma above we can have an idea how the graph of g{r) looks like: 




Define Kq = y/ g{rQ) = yf^-^ ^('^o)- From the Lemma follows that the restriction of g{r) 
to the interval (IjJ'o] is a decreasing function. Also g{r) restricted to the interval [ro,oo) is an 
increasing function. Hence the same is true for the function y/g{r). Therefore the functions 
y/g\(i^ro] '■ {^ii"o] [KqjOo) and y^|[ro,oo) ■ [''0,00) — > [Kq,oo) have inverses which we call ri 
and r2 respectively. It is not difficult to verify that {iv), {v), (vii), (viii), {ix) hold for ri and 
r2. Next we prove {vi). Fix K > Kq. Note that U{r) = ^ + V{r) ^ ^ ~ 4A/q^ ^^^Ztn^e ^ 

- Since - 4A/o^^^i — ) tends to -27rA when r ^ oo, then for 

large r we have that U{r) is negative. This together with (ii), {Hi) of the statement of Theorem 
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C and the fact that U{r) has exactly two critical points at ri and r2 imply that ^ U{r) > for 
1 < r < ri and r > r2. To prove that ■^U{r) < for ri < r < r2 it is enough to prove that 

A U{ro) < 0. But i U{ro) = + i V{ro) = -^^ + ^ = < 0. 

'o 'o 'o 'o 

Finally we prove (x). Let f denote ri or r2. Since /C^ = r^^{f) we have ?7(f) = ^ + 

2-^ /o^ (f2!feof2^)3/2 ^^i where the last equality is obtained by a change of variable. 

If — > +00, we have (by {viii) and (ix)) that ri and r2 — > +oo. Hence it is enough to 

prove that lim^^i /^^ ( ^^"^^^fa^'p/a dO = +oo and lim^^+oo Jo'' j^^^^^dO = 0. The last hmit 
is clearly zero. We prove now Muir^i ^ (.J^f ^2^)3/2 dO = +00. We have J^^^^^dO = 
/o' ( j!_°cof2e)3/2 + /f (r2-co!2e)3/2'^^- hm^^^i /| j^^^^^^dO = /| "^""-Jg^dO is finite, 

it is enough to prove that lim,._+i Jq'' ^^2Tcog2^)3/2 = 00. We can assume r < y^. Hence 
2cos^e-r^ > i, for 9 G [0, f ]. Therefore /| > 1 / 1 . But 

since (r^ — 1) < |, sin^O < sinO < 1 and sinO < 9 we have (r^ — 1) + sin^O < |. Hence 
((r2 - 1) + sin^ef^ < ^((r2 - 1) + sin9) < ^((r^ - 1) + 0). 

Consequently ^ ((,2_i)+l2g)3/2 > ^ /o^ = ^[/^((r^ - 1) + f ) - ln{r^ - 1)] ^ 00 

as r — > 1. ■ 



Before we prove Proposition 3 we need some definitions and comments. Since V{r), r G 
IR^ — C, is invariant by rotations around the z-axis we can reduce our problem in a canonical way 
to a problem with two degrees of freedom. Using cylindrical coordinates (r, (p, z), the Lagrangian 
in these coordinates can be written as: L{r,ip, z,r,ip, z) = ^(r^ + r'^ip'^ + i^) — V{r,z), where 
V{r,z) = V{t), r = (r cosip,r sirup, z). Then it is straightforward to verify that the system 
r= — Vy(r) in these coordinates is given by: 



— 



dV 
dr 



dv 

az 
K 



(5.7) 



where K is the (constant) angular momentum. Note that the first two equations of system (j5.7j) 
can be rewritten as {r,z) = —VU{r,z^, with u{r,z^ = ^ + V{r,0,z). If {r{t),z{t)) is a 

solution of the first two equations of (|5.7|) . defining ip{t) = we have that (r(t), p{t), z{t)) 

is a solution of 1)5. 7p . Then {r{t) cosp{t),r{t) sinp{t), z{t)) is a solution of r^= — Vy(r). Note 
that if z = 0, system (|5.7j) is reduced to (|3.4|) and [/ becomes the effective potential U. 



Proof of Proposition 3. Suppose s 
Vl7(s) = 0. In particular = 0. But |^ = where (3 

solution lies in the horizontal plane. Since z 
effective potential. 



{r,z) is an equilibrium position of s= — V{7(s). Then 
Jc ||r_tj|3 > 0. Hence z = 0, i.e the 
C/(r), where U is the 



C lr-u\\^ 
0, we have that U (r, 0) 
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Let (r, 0) be an equilibrium position of s= — Vt/(s). If 1 < r < ro we know that the 
corresponding circular solution is not stable. To prove that the circular solution that corresponds 
to (r, 0), r > ro, is stable it is enough to prove that (r, 0) is a strict local minimum of U. 

Since (r, 0) is a critical point of U, we have that r is a critical point of U, and since r > tq, 
r = r2(|-R'|) where K is the angular momentum of the solution. Hence r is a strict local minimum 
of U (r) that is U{r) < U (r') for r' close to r, r' ^ r. But an easy calculation from the definition 
of V shows that V{r',0,0) < V{r',0,z), for all r',z, {r',z) ^ (±1,0) and the same holds for U. 
This implies that U{r', z) > U{r', 0) = U{r') > U{r) = U{r, 0) for r' close to r, r ^ r' . Therefore 
r = r2(|-fC|) is a strict local minimum ofU . m 

To finish this section we prove that the origin is the only equilibrium position. Note that 
by symmetry it is a simple matter to show that the origin is in fact an equilibrium position. 

Proposition 5.2 The origin is only equilibrium solution of the system r= — W{r). 

Proof. For r = {x,y,z) differentiating we have that ^(r) = Pz, where P = Jc jr-up 
Then if W{r) = then z = 0. Hence every equilibrium position lies in the a;y-plane. The 
Proposition now follows from (i) of Theorem B and {i) of Theorem C. ■ 



6 Proof of Proposition 4. 

Fix A > 0. Consider the fixed homogeneous circle with constant density A, contained in 
the xy-plane and passing through the origin, with radius ^ and center (^,0,0). Note that the 
xz-plane is an invariant subspace of our problem. Denote the potential of this translated fixed 
homogeneous circle, and restricted to the xz-plane, by W{x,z;e). For notational purposes in 
what follows we use coordinates (x,y) instead of coordinates {x,z). 

This potential can be written in the form W {x,y;e) = — 4A-T(D,d) where T,D,d are as 
in section 1.1. If x < we have that = (f — x)"^ + v"^ and cP = x"^ + y"^. Hence we have (see 
section 1.1) W (x, y; e) = -4A§/ (^j^^^j^) /^(c-x)^ + where c = f . 

Let "VW denote the gradient (^, ^) of W. A straightforward calculation shows that for 
X < c = I we have 

Lemma 6.1 

VW (x, y; e) = 2Ac^(x - c,y) - 4A-^/'(0 [c(x, y) + (y^ - -2xy)], 
where t = ^^llyily-i and D = {c - x)'^ + y'^ . 
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Remark. If z = x + then = (x^ — 2x?/), hence the identity above becomes 
VW {z; I) = 2Xj^f ( J^) (z-c)- 4A^/' ( J^) {cz - z') = 

= -2A^{/ (^) + z (^) /' (^) } (. - c). 

Note that the gradient vector has components in the z and z^ directions and in the real (1,0) 
direction. 

, / N 4A/(t) 1 , , , 32A 1 1 , , , , 16A/'(t) 1 
Now set hi{x,y-e) = h2{x,y]e) = - — j^-^ and h3{x,y;e) = — 

Hence 

VW{x,y;e) = [ehi]{x,y) - [2^ehi]{l,0) + [e^f'{t)h2]{x,y) + [Ve hs]{x^ - y\ 2xy) . 
The next Lemma shows that hi,h2, /13 are bounded in certain sense. 

Lemma 6.2 Given a compact C C IB? — {0} there are eo > 0, X > such that \hi{x, y; e)| < K 
for all (x, y) € C, < e < eo- 

Proof. Given a compact C C IR^ — {0}, there exist constants Ci,C2, such that < Ci < 
||(x,y)|| = \/x2 + 1/2 < (^2, for all (x,y) G C. For eo < we have C2 < which implies that 

C C {(x,j/);x < -^,2/ € M}, for all e < eo, and hi{x,y;e) is defined in C, for e < eo, i = 1,2,3. 

Since V(2-ex)2 + eV = ||(2, 0) - e (x, y)|| < ||(2, 0)|| + e || (x, y)|| = 2 + e\\{x,y)\\ and 
11(2,0) -e(x,y) II > ||(2,0)|| - e ||(x,y)|| = 2 - e ||(x, y)||, and we obtain 



2-eC72 < ^(2-ex)2 + eV < 2 + eC2 (6.8) 
Therefore, using Lemma ll. 21 part (iii), and (|6.8() we have that: 



\hi(x,y-e)\ = 



4A/(t) 1 

e D3 £3/2 



((2-es)2 + e2j;2)3/2 - (2-eC2)3 " (2 - eCi)^ ^ 



^ 2nX^[{2-exf + e%Y^^ < 27rA(2 + 6^2)^/' ^ 2^A(2 + = ^ 

(2-eC2)«,/^(:r2 + y2)i/4 - ^2-^02^01^' " (2 - eoC2)3<' 

This proves the Lemma for hi. 

Now, applying (jSSl) we have \h2{x, e)| = = ((2-..)'+e2y^)V2 < jT^lkv ^ {2-clk}-' = ^2(C, eo), 
and this proves the Lemma for /i2. 

Finally, using Lemma lTT^ part (v) and ()6.8() . we have |/i3(x, y; e) 



16A/'W 1 


16A 


l/'WI 






— ^5/2 


(2-ea;)2 + e2j,2 


5/2 






.2 





_ 16Ae5/2|/'(f)| 4A7r 1 ^ 4\n 1 ^ 4A7r 1 _ i^- , N _ 

((2-ex)2+e2y2)5/2 - (3,2+^2)5/4 ((2-ez)2+e2y2)5/4 ^ (2-eC2)^/2 - cf^^ (2-eoC2)5/2 ^3{'^,t0)- ■ 

Write /i(x,y;e) = {x'' + y^)e^f'{t)h2. Then VT^ = [ei/ii](x, y) - [2^^/11] (1, 0) + /i (j^ij^ + 
[Ve/i3](2;^ -y^2xy). 
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Lemma 6.3 // {xn,yn,(-n) {x,y,0) / (0,0,0), with e„ / 0, (x„,y„) ^ Ce„, then 
limn^+ooh{xn,yn;en) = 64A. 

Proof. First note that lim„^+oo h2{xn,yn;^n) = -32A. Also lim„^+oo = -^iip^, where t„ = 

- — ^^'^ — -. Therefore hm„^+oo /i(a;„,y„,e„) = hm„^+oo(a:^+y^)(7^)(/'(in)in)/i2(a;n, yn; Cn) = 
(^"-^j +yl 

(x^ + ^)(— 32A) = 64A, where we are using the fact that \mii__,Q+ t f {t) = — ^ (see 

Lemma ll.2p . ■ 



Define VW{x,y;Q) := 64A^#i^, (x,y) / (0,0). Hence VW(x,y;e) is defined on ^ = 

X -\-y 

{ (x, y; e) ; + 7^ 0, (x, y) 7^ (7,0) }, and for each e, VW is analytic. 

Proof of Proposition 4. Clearly VW is continuous in (x, y, e), with e 7^ 0. Let { {xn,yn, ^n) } 
be a sequence in A, with (x„, y„, e„) ^ (x, y, 0), (x, y) / (0, 0), e„ 7^ 0. Thus C = { (x^, y^), n E 
IN } U { (x, y) } is compact and (0, 0) ^ C By Lemma 16.21 there exist no and Kq > 0, such 
that |/ii((x„,y„,e„)| < Kc, for n > no, i = 1,2,3. Therefore, lim„_^+oo VW^(xn, y^, e^) = 
lim„_+oo/i(x„,y„;e„)5?^ = 64A-fe4. ■ 



7 Proof of Proposition 5. 

Take M = 1. Let /3(t) = {D{t),d{t)) with /3(0) = (Z:',d), /3(0) = {u,v). Let us calculate 
o"(t) = a{D(t),d{t)). The following change of variables will simplify our calculations. Write 
u = and w = (7 + 5)d. The Taylor series of D{t), d{t) around t = are: 

D{t) = D + tu + ... = {1 + t-f)D + 0(^2) 
d(t) = d + tt> + . . . = (1 + t7)d + t5d + O(t^) 



Lemma 7.1 f7(t) = (l+ta*)fT + 0(t2), with a* = ao + 5J2t=i ^^U]Zl (l 
do = d, Do = D, Dn = ^"-iW+'^-iW ^ dn = V^„_i(0)d„_i(0), /or n > L 

Proof. Calculating the arithmetic-geometric mean of D{t) and d{t), we have: 

^ D{t) + d(t) ^ ^ ^^^^^ ^ 1^^^ ^ ^^^2) ^ ^ ^^^^^^ ^ ^^^2)^ 



; "0 = 7; 



di{t) = \lD{t)d{t) = (1 + t7)di + -tSdi + 0(r) = (1 + tai)di + -tWi ( 1 - ) + 0(r 



where ai = 7 + i^;^, = ^^2)+m. = D±d ^ /D(0)d(0) = 

Calculating the arithmetic and geometric mean successively, rearranging and defining 
ao = 7; ai = 7 + and 



23 



we obtain, by induction, the general rule: 

dn{t) = (1 + tan)dn + —tSd^ [] (^1 " ^ j + 0{t^), 

where do = rf, Dq = D,Dn = W+'^-i W , d„ = ^Dn-mdn^iW, n>l. 
Note that an = ao + d ELi n}=i (l - • 

Finally, taking the limit on dn{t) or on Dn{t), when n goes to infinity, we have 

a{t) = {l+ta*)a + 0{t'^), 

with a* = lim„_>oo an = ao + 6 J2t=i 2^%^ ^l=i ~ "TT^) • ^^^^ ^^^^ series 

converges. ■ 



Proposition 5. We have the following formulas: 

^V{D, d) = ^V{D, d), ^V{D, d) = -^V{D, d), 

where x = E^n=i 2^%^ UU (l " %f ) • 

Proof. Since V{P{t)) = ^ = i - *^ + 0{t^) = V{D,d) - ta*V{D,d) + 0{t^) we can now 
calculate the partial derivaties -^V{D,d), -^V{D,d). 

We have j^V{/3{t))\t=o = -a*V{D,d) = ( (^y(/3(0)), ^F(/?(0))), {u,v)). Note that 
a* = 7+<5x, with X = E+=i 2^%^ YT.Zl (1 - %f ), 7 = ^ and <^ = 1-^. We have a* = § + 
{2-%)x = ^(l-x)+ix. Therefore a* = ( V) so ( (^Fp, d), ^V{D, d)), (u, v) ) = 

- (^^V{D,d)u+ lV{D,d)v) for ah u,v. Hence: 

= {^v{D,di m^)) •■ 
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